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Abstract. We present a framework for studying gravitational lensing in 
spherically symmetric spacetimes using 1+1+2 covariant methods. A general 
formula for the deflection angle is derived and we show how this can be used to 
recover the standard result for the Schwarzschild spacetime. 



1. Introduction 

In this paper, we apply the recently developed 1+1+2 covariant approach to the lensing 
of light rays in any spherically symmetric spacetime. A more detailed discussion of 
this approach can be found in [5], from which most of the material in section [l"T1 has 
been drawn. We begin with a motivation for using the 1+1+2 approach rather than 
the conventional 1+3 covariant method used in a related paper by the same authors 
to derive the cosmological lens equation pQ. 

The 1+3 covariant approach [2] [3] has proven to be a very powerful framework in 
cosmology, not only in studying the motion of null rays and gravitational lensing, but 
in many other areas of cosmology. In particular this approach has been particularly 
useful in obtaining a deep understanding of many aspects of relativistic fluid flows, 
particularly in areas involving perturbation methods. In cosmology these methods 
have been applied to the evolution of linear and non-linear fluctuations in the universe 
and to the physics of the cosmic microwave background [4], amongst other things. The 
strength of this approach lies in the fact that all the essential information describing a 
particular system is captured in a set of 1+3 covariant variables which are defined 
relative to a preferred time-like congruence, often taken to be the 4-velocity u a 
of matter, allowing these variables to have an immediate physical and geometrical 
significance. These variables then satisfy a set of evolution and constraint equations 
which are derived from the Einstein Field Equations, the Bianchi and Ricci identities, 
forming a closed system of equations for a chosen equation of state describing the 
matter. 

A natural extension of the 1+3 approach involves a further splitting of the 3-space. 
This 1+1+2 decomposition of spacetime is ideal for solving problems which have 
spherical symmetry. In particular this approach was used to study linear perturbations 
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of a Schwarzschild spacetime and to investigate the generation of electromagnetic 
radiation by gravitational waves interacting with a strong magnetic field in the vicinity 
of a vibrating Schwarzschild black hole [5] . 

In addition to the splitting with respect to a timelike vector u a , this 1+1+2 
approach relies on the further decomposition of spacetime using a spacelike vector e a . 
The Bianchi and Ricci identities can then be split using both u a and e a giving rise to a 
set of coupled first order differential equations and constraints. Differential operators 
are along these two vector fields resulting in a set of evolution and propagation (along 
e a ) equations (see [5] for the full set of equations). We also denote a differential 
operator orthogonal to the spacelike vector e a (and thus lying in the 2-dimensional 
sheet) by S a . 

The 1+1+2 decomposition of the null tangent vector k a gives more insight on the 
lensing geometry for a light ray traveling in a spherically symmetric spacetime and 
allows us to directly obtain the general form of the deflection angle a, for spherically 
symmetric lensing situations, which is the main result of this paper. 

Conventions are taken to be the same as in [1] and unless otherwise stated, we 
use geometrized units, so that 8nG = c = 1. 

1.1. The 1 + 1 + 2 Covariant Sheet Approach 

In the 1+3 covariant approach, a 4- velocity field, u a = dx a /dr, is introduced which 
represents the average 4- velocity of a set of observers in the spacetime. The projection 
tensor, h ab — g ab + u a Ub, can be used to project any vector or tensor orthogonal to 
u a . Thus, h ab allows for the decomposition of any 4- vector into a component parallel 
to u a and a component lying in the 3-space orthogonal to u a . 

Another vector field can be defined which allows us to perform an additional 
decomposition of the 1+3 equations. Introducing a unit spatial vector e a which lies 
orthogonal to the u a so that u a e a — and e a e a — 1. The projection tensor 

N a b = h a b - e a e b = g a b + u a u b - e a e b , N\ = 2 , (1) 

then projects vectors orthogonal to e a and the timelike vector u°, i.e. e a N ab = = 
u a N ab . Thus, N ab projects a tensorial object onto 2-surfaces which we call the 'sheet'. 

In the 1+3 covariant approach to the study of null congruences [Tj, we introduced 
the spatial vector, n a , which was chosen such that n a n a = 1 and u a n a — (i.e. n a 
is a unit spatial vector). The tensor, h a b = h a b — n a nb, could then be defined as 
a projection tensor, which projects a vector or tensor into the 2-dimensional screen- 
space, orthogonal to the null tangent vector k a . In the following analysis, we have 
not allowed for the spatial vectors n a and e a to coincide. However, if we choose n a to 
coincide with e a , it immediately follows that h a b — N a i>, so that the screen and sheet 
represent the same 2-dimensional surface. 

Any 3- vector, ^ a ', can be irreducibly split into a component along e° and a sheet 
component orthogonal to e° i.e. 

tjj a = q> e a + * a , * = iP a e a &nd * a = N ab ip b . (2) 

A similar decomposition can be done for a projected, symmetric, trace-free (PSTF) 
tensor, ip a b, which can be split into scalar, vector and tensor part as follows: 

^ab = i>(ab) = *(e a e b - ^N ab ) + 2* (a e b) + * ab , (3) 
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and the curly brackets denote the PSTF part of a tensor with respect to e a . We also 
have 

2 

h{ab} = , N {ab} = -e< a e 6 ) = N ab - -h ab . (5) 

We now define the alternating Levi-Civita 2-tensor as 

£ab = e abc e c = rj dabc e c u d , (6) 

where e abc is the 3-space permutation symbol and rj abc d is just the spacetime 
permutator. With the definition of e ab given by equation ^ above, we also have 
the following relations : 

£abe b = = £( Qb ) , (7) 
£abc — ^a^bc ^b^ac H~~ &c&ab ? (8) 

e ab e cd = N a c N b d - N a d N b c , (9) 

e a c e bc = N ab , e ab e ab = 2 . (10) 

From these definitions it follows that any object in the 1+1+2 setting can be split 
into scalars, 2-vectors in the sheet, and PSTF 2-tensors (also defined in the sheet). 
These are the three objects that remain after the splitting has occurred. In the 1+1+2 
formalism we can introduce two new derivatives: 

^ a .. b c - d =efV f ^ a .. b c - d , 

6fip a .. b c ~ d = Nj ...N b g N i c ...N j d NfV j i> f .. g i ^ . (11) 

The hat-derivative represents the derivative along the e a vector-field. The S -derivative 
is the projected derivative onto the sheei[§|. 

We are now able to decompose the covariant derivative of e a orthogonal to u a 
giving 

D Q eb = e a a b + -(j)N a b + &ab + Cafe > ( 12 ) 

where 

a a = e c D c e a = e a , (13) 
4> = S a e a , (14) 

£ = ±e ab S a e b , (15) 

Cab = S{ a e b } . (16) 

| Note that for any 2-vector *f a , we may use £ a j, to form a vector that will lie orthogonal to \E" a but 
of the same length. 

§ We note that one needs to project on every free index when calculating the 5-derivative. 
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We see that when travelling along e a , <f> represents the expansion of the sheet, £„& 
is the shear of e a (i.e. the distortion of the sheet) and a a its acceleration. We can 
also interpret £ as the vorticity associated with e a so that it is a representation of the 
'twisting' or rotation of the sheet. The full covariant derivative of e a is 

V a e 6 = - Au a u b - u a a b + ( E + -6 ) e a u b + [£„ - e ac il c ] u b 



+ e a a b + -(j>N ab + £e ab + ( ab 



(17) 



and 



V a M b 



u a (Ae b + A b 



e a e b 



-e 



(Z a - e ac n c ) e b + N ab [ ie-^E 



e a (E 6 + £ft c fi c 



fie a ft + E a ft , 



(18) 



for the decomposition of the 1+3 covariant derivative of u a . These derivatives will be 
used in later calculations. The derivative of e° in the direction of u a is given by 

e a = Au a + a a , where A — e a it a , (19) 

and a a is the component lying in the sheet. The new variables a a , <fi, £, £ a b and 
a a are fundamental objects of the spacetime and their dynamics give us information 
about the spacetime geometry. Essentially, they are treated on the same footing as 
the kinematical variables of u a in the 1+3 approach. 
The covariant derivative of a scalar \1/ is 

D o tf = $e« + * tf. (20) 

While for any vector ^ a orthogonal to both u a and e° (i.e. ^ a lies in the sheet), the 
various parts of its spatial derivative may be decomposed as follows [jj : 



D Q * b = -e a e b V c a c + e a ^ b 



^a + {S,£ac + Cac)^ C 



Similarly, for a tensor "i> ab (where ^ a 6 = ^{ab}) 

"1 



D a *6c = -2e a e {b ^ c)d a a + e a * bc - 2e 



^c)a + *c) d (&ad + Cad) 



be 



(21) 



(22) 



We include the following relations which may be useful: 

Nab = 2w( Q u 6 ) - 2e (a e fc ) = 2u( a A b ) - 2e( a a b ) 



N ab = - 
S c N ab = 



2e( a a b ) 



(23) 



while 



e ab = — 2u[ a e b ] c A c + 2e[ Q e 6 ] c a c , 
iab = 2e[ a e b ] c a c , 

ScSab = . 

Note that a bar on a particular index indicates that the vector or tensor lies in the sheet. 



(24) 
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We now split the (1+3) kinematical variables and Weyl tensors into the irreducible 
1+1+2 set {6, A, O, E, S, H, A a ,Z a ,£ a ,H a , E q6 , £ ab , H ab } using equations © and © 
giving: 

ii a = Ae a + A a , (25) 
uj a = fie a + f2 a , (26) 

CT a b = E (e a e b - 7j N ab^J + 2S (a e b) + E ab , (27) 
^6 = £ [e a e b - \N a ^ + 2£ {a e b} + £ at> , (28) 

H ab = H [e a e b - ^N ab ^j + 2U (a e b) + U ab . (29) 

We have now obtained the 1+1+2 form of the kinematical and Weyl variables 
corresponding to the irreducible parts of V a e (|17j) . Similarly, we may split the 
anisotropic fluid variables q a and TT ab '- 

q a = Qe a + Q a , (30) 
TTab = II (e a e b - 7; N ab^J + 2iT( Q e fc) + IT a & . (31) 



2. Gravitational Lensing in Spherically Symmetric Spacetimes 

This section begins with a presentation of the various kinematical relations valid in 
the general local rotationally symmetric (LRS) spacetime. A detailed discussion of 
the covariant approach to LRS perfect fluid spacetimes can be seen in [7]. However, 
in deriving the lensing geometry we have restricted ourselves to the case of spherically 
symmetric spacetimes. We also derive the propagation equations for the lensing 
variables along k a in the general LRS spacetime. These equations will be needed 
in section [4] where we derive the deflection angle in the Schwarzschild spacetime. 



2.1. Kinematical Quantities 

As discussed in sectionQ] the 1+1+2 approach requires a further splitting of spacetime 
variables, in addition to the split with respect to the time-like vector u a . This extra 
splitting was done relative to the vector e a . In the case of LRS spacetimes, e a is just 
a vector pointing along the axis of symmetry and can therefore be thought of as a 
'radial' vector. LRS spacetimes are defined to be isotropic, allowing for the vanishing 
of all 1+1+2 vectors and tensors, so that there are no preferred directions in the sheet. 
The full covariant derivative of the radial vector e a (fT7| and 4-velocity u a (fTBj) reduce 
to 

V a e h = -Au a u b + ( S + ^9 ) e a u b + \<t)N ab + £e ab , (32) 



and 



3 



V a M 6 = ~Au a e b + e a e b ( ^6 + S ) + N ab ( ^9 - I + fle ab , (33) 



3 J V 3 2 
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respectively. We also find 

e a = Au a , (34) 

is the evolution of e a . 

The kinematical quantities and Weyl tensors ([25]) -([3Tj) reduce to 



u a = Ae a , 


(35) 


cu a = ne a , 


(36) 


1 

<7ab = E(e a e b - -iVafc) , 


(37) 


q a = Qe a , 


(38) 


Tab = II (e a e b - ^N a bj , 


(39) 


E a b = f(e Q e fc - -N a b) , 


(40) 


Hab = H(e a e b - -Nab) , 


(41) 



in an LRS spacetime where we have not made the perfect fluid assumption (i.e. 

O^n^o). 

We now list the set of 1+1+2 equations resulting in any LRS (and thus any 
spherically symmetric) spacetime: 

\&-e)(a-U) +2& + Q , (42) 



-v--Q)t + \A--*)n + -H, (43) 

n = +At + n (z- ^ej , (44) 

A = ^9 - £ - A 2 + Qe + + i (/i + 3p - 2A) + £ - in , (45) 

i = - (u + 0) .4 + ie 2 + ^£ 2 - 2Q 2 + ~ (p + 3p) - A , (46) 

2-1 /2 1 \ 2 „ 1 

£ = -A + - (2A - </>) u - I -9 + -£j £ - -O 2 - £ + -IT , (47) 

A = - Q - 6 (/x + p) - (0 + 2A) Q - ^£11 , (48) 

q= -p-n- f^+^Jn- Qe + EjQ-(^+p)^, (49) 



2 3 V 2 J 2 V 3 
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i(i0-2u)Q + 3£«-i(/i+p)E, (50) 



H= -HE+ 6 + -S )H + CLQ + -£Il , (51) 



1,2 /1„ „\ /2 



- ~(/x + A)-~n-£, (52) 
£ = -#+(~e + s)n, (53) 

E - ^© - - - Q , (54) 
CI = + (A - <£) fi , (55) 

^ = + + ^ n ) +l^Q + 3nn > ( 56 ) 

?? = - \<\>n- Us + v + p- hij n-Qt . (57) 

= (2A -(f>) CI — 3££ + % (58) 
2.2. Lensing Geometry 

In pQ we introduced the spatial vector, n a , which is a 3-vector in the direction of light 
propagation. Allowing for a further split of this vector with respect to e a gives : 

n a = ne a + n a . (59) 

Thus, k is the magnitude of the radial component and k° is the component lying in 
the 2-dimensional sheet. In this way we are able to write the null tangent vector k a 
as 

k a = E(u a + ne a + n a ) , (60) 

where E = —k a u a is the energy of a photon relative to u a . Using this form of the null 
vector, we can determine the lensing geometry of a photon experiencing a deflection 
about the center of symmetry. 



2.3. General Form of the Deflection Angle 

Figure la traces the path of a photon in the presence of a strong gravitational field 
such as a black hole, with k a lying tangent to the null geodesic. We can also see 
the directions of the various vector components of the null vector. The 2-dimensional 
sheet is represented as a sphere about a central point. When considering infinitesimal 
small angles of deflection da, the lensing geometry will be given by figure lb. The 
angle da subtends an infinitesimal small displacement of the photon's path, given by 
|i?K a |. da can be obtained from figure lb giving 

da = - \En a \ dv , (61) 
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Figure 1. Lensing geometry in symmetric spacetimes 



which may be integrated giving the total deflection angle as 

a = I - \En a \dv-a , (62) 

J u± ' 

where ao is the integration constant. This is a completely general form of the scalar 
deflection angle which can be applied to any spherically symmetric spacetime. The 
deflection angle a ([62]) may be written as 

r"2 



a 



a 



f" 2 1 , 

/ - \E\ ^lT a dv - a , (63) 

J u 1 ' 

where the magnitude of n a (i.e. \n a \) is just ^k"k„. Substituting for k a (foP]) into the 
null condition, k a k a = 0, gives a restriction on K a : 

n aK a = 1 - k 2 , (64) 

so that the general deflection angle takes the form: 

/ - |J5| Vl - &du - a Q . (65) 

J V\ ^ 

Thus, if we know the relations E(y), k(v) and r(u), for a given spherically symmetric 
spacetime, it is possible to find an explicit form of the deflection angle. The differential 
equations for these variables in the direction of the ray are derived in the next section. 

2.4- General Propagation Equations along k a 

We begin this section by deriving the propagation equations for the lensing variables 
E and k in the direction of the null ray. Each ray is parameterized by the affine 
parameter, v 1 so that the geodesic condition can be written as : 

fik a 

k b V b k a = ?L = (k a )' = , (66) 
dv 

where the prime derivative (') denotes change along the ray (i.e. with respect to 
the affine parameter v). The geodesic condition (|66)) can be used to determine the 
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propagation equations for E and k. Substituting for the null vector k a (f6"0"|) into (f6"6")) 
gives 

1 









U a + 





E' + -E 2 <h + E 2 (-0 - -S 
2 V3 2 



where we have also used 



1 



£"«; + £ 2 „4 + £ 2 k [ -9 + £ ) + Ek' 



K a + E 2 (n + e Qb K b + = , 



1„ 1. 



(67) 



< = EAe a + Ek -6 + £ e + £ -6 - -£ ) K a + Ene ab n 



e' a = EAu a + Ek ( E + -6 ) u a + --E^«; a + E£e ab K b , 



(69) 

which are obtained from (fT8|) and (fT7|) . respectively. When (|67|) is projected along the 
timclike direction (u a ) and radial direction (e a ), we obtain the general equations of 
E' and k', respectively. These are : 



E' = - E 2 kA - 2E 2 k 2 T i - E 2 ( -O - S ) , 



k' = - £7 (1 - k 2 ) \ A- -(/>+ 2kS ) , 
which can be used in a general LRS spacetime. 



(70) 
(71) 



The LRS class II spacetime [7], admits spherically symmetric solutions. This space- 
time is free of rotation allowing for the vanishing of the variables H, fi and £. If in 
addition we describe the matter to be a perfect fluid so that Q = II = 0, then we can 
determine the propagation of some of the scalars, {/i, p, 0, S, A, £, </>}, along k a 
usin^ 

X' = E (X + kX\ , (72) 

where X can be found through the addition and subtraction of the relevant 1+1+2 
equations (|^T) - (t55|) Fl. 



3. Application to Lensing in Schwarzschild Spacetimes 

Schwarzschild gravitational lensing in a weak gravitational field (i.e. where the 
deflection angle is small) is well known (see for example [5]). An extension of 
Schwarzschild lensing in the strong field limit is investigated in [TU] and [5]. 

We have general forms for the deflection angle and the propagation equations 
of the lensing variables along k a . These will now be applied to the Schwarzschild 
spacetime in which we hope to recover the standard lensing results found in [8]. 

1 Note that 5 a X = in a spherically symmetric spacetime. 

We also note that in using a spherically symmetric spacetime, these equations will reduce to simpler 
expressions as"H = f2 = £ = and Q = II = 0. 
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3.1. Solutions for A, £ and <f> 

The Schwarzschild spacetime is characterized as being static, spherically symmetric 
and a vacuum. This covariant characterization of the spacetime allows for the 
quantities A, £ and <f> to be non-zero, with all the other covariant quantities vanishing. 
The remaining propagation equations are [5] 

j>=-^cb 2 -£, (73) 

i=-^4£, (74) 

together with an extra equation 

£ + Ac/) = . (75) 

These are a closed set of equations from which we are able to find solutions for the 
variables A, <f> an d £■ The coupled differential equations seen in (|73l) and ([74]) can be 
solved explicitly when the hat-derivative is associated with an affine parameter, say 
p, so that~= d/dp. The parametric solutions are : 

sech 6 x = — , (76) 



(2m) 



1 2 / 2m 

4> = — sink x sech 3 x = —\1 , (77) 

m r V T 

1 , m / 2m \ 1//2 

A = - — cosech x seen x = — 1 , (78) 

4m r \ r J 

where solutions can be written in terms of either parameter x or r, and are related to 
the affine parameter through the relation 

p = 2m (x + sinh x cosh x) . (79) 

r is the usual Schwarzschild radial coordinate related to parameter x by 

r = 2m cosh 2 x . (80) 

Solutions (1751 - ([75)1 are a one-parameter family of solutions, parameterized by the 
Schwarzschild mass m. Note that the Schwarzschild solution is given for 2m < r < oo 
for < x < oo. 



3.2. Solutions for the Lensing Variables 

As stated in the previous section, the only non-zero 1+1+2 variables in a Schwarzschild 
spacetime are the scalars {A, £, <fi} (and their derivatives {A, £, 4>}). Thus, the general 
propagation equations of E (|70"|) and k (I7T|) . in the direction of the ray, reduce to the 
form : 

E' = - E 2 Ak , (81) 
k' =E(l- K 2 )(^-A) . (82) 

* This covariant characterization of the 1+1+2 variables in a Schwarzschild spacetime is also 
discussed in [5]. 
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Equation (|64[) gives the propagation equation for k" when differentiated with respect 
to v : 

K al = -En \-<j> - A ) n a - -E<j> (1 - n 2 ) e a . (83) 



In order to find solutions to these differential equations, we also include the 
propagation equations for A, £ and (j). These are obtained by substituting for the 
1+1+2 form of k a (|60|) into the definition of the prime derivative of a particular 
variable. We then find: 

A' = k a V a A= -Ek{(/) + A)A , (84) 
0' = k a V a( j) = En(A - ^cj>)<j> , (85) 

£' = k a V a £ = ~E K <f>£ . (86) 



The differential equations (|M]> -(|86 | can be used in solving for E and k as will be 
described below. 

Dividing equation (1551) by </> gives: 

— = EkA - \ek4> . (87) 
(f> 2 

We can then use the differential equations for E (|8Tj) and £ (|86|) to rewrite (f87| as: 

E' _ ft 1 £' 

~E ~ + 3~£ ' 
which can be solved as a total differential equation, i.e. 



(In E)' = (in <f> - In £ 1/3 ) = ( In -| 



£1/3 



=^3, (89) 

with Ci being an integration constant with respect to affine parameter v. Similarly, 
we are able to solve for k starting with (|85j) and using the differential for k (|82|) . We 
find: 

<b' , . 

5 = ~ , (90) 



which can be written in total differential form as 
1 



t\n (1 - k 2 ) 



= (\n <f>) ' . (91) 



This equation can be integrated giving 

>-£) 4 : 

where Ca is a constant of integration with respect to v. 

We can now rewrite the solutions for E and k in terms of the radial coordinate 
r and Schwarzschild mass m. In order to do this, the integration constants C\ and 
Ci need to be determined. C\ can be obtained by taking the limit of equation 
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as r tends to oo, and given the solutions of £ and <p from equations ([75)) and (|77|) 
respectively. We find 

(2m) 1 /3 



Ci 



2£ D 



(93) 



where we denote the energy at a large radial distance as, E^. 

At the point of closest approach 0, r — r , we must have k — so that we are 
able to determine Ci from equation ([9"2")) above. We obtain : 



C, = 



/4(r - 2m) 



(94) 



where ro is the distance of closest approach. Substituting for the constants of 
integration ([9"3")) and (|M|) . into ([59")) and ([9"2"]l respectively, gives : 



E = E n 



1 



2 m 
r 



(95) 



± 



1 



1 



(96) 



r 2 (r — 2m) 

Equations ([9"5"]) and together with the solutions for {.4, £ </>}, form a complete 

set of solutions from which the lensing geometry in a Schwarzschild spacetime can be 
determined. This is described in the next section. 



4. The Deflection Angle a 



We would like to rewrite the scalar deflection angle ([6"5")) in terms of the Schwarzschild 
mass m and radial coordinate r so that it may be compared to that given in the 
standard lensing literature. Substituting for the solutions for E ([95]) and k ([96]) into 
(f65|) gives 

1/2 



U 2 E 

r 



r 2 (ro — 2m) 



dv — otQ 



l '2 



J 



where J is the impact parameter defined as 



2m 

J = r 1 

7-0 



-1/2 



(97) 



(98) 



We need a transformation relation between the affine parameter dv and radial distance 
dr. Differentiating the solution for £ ([76"]) with respect to v gives: 

6m , 



£' 

and substituting for k ([96 
results in 

6m 



(99) 



77)) and £ ()76l) into the differential equation for £ 



2m 



1/2 r 



1 - 



r 2 (ro — 2m) 



1 - 



2m 



1/2 



(100) 



(I Here rg is just the closest point that the ray would reach in the vicinity of the lensing object. 
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The transformation relation can be obtained in equating equations (IM|) and (|100D 
giving 



dr = E a 



1 



ro" 



1 - 



2m 



E a 



r 2 (tq — 2m) \ 

1/2 



1/2 



dv , 



fi/./ 



(101) 



2m 

Using this transformation in equation (|97|) , gives the form of the scalar deflection angle 
in terms of the impact parameter J and radial coordinate r as : 

-1/2 

dr - a , (102) 



a = 2 



J 



J 2 / 2m 
1- — [1- 



which is the exact expression obtained in the lensing literature [8] where a metric 
based approach is used (and if we take ao = tt)- 



5. Conclusions 

In this paper we have developed a new approach for studying lensing effects in 
spherically symmetric spacetimes based on the 1+1+2 covariant approach to general 
relativistic fluid dynamics. In particular, we derived a completely general formula for 
the bending angle which is expressed in terms of the geometry of the null congruence. 
This approach proves particularly useful in the study of lensing in f(R) theories of 
gravity QT]. 
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